N73 24936 

1 

1 

THERMODYNAMIC PROPERTIES OF UF g 
MEASURED WITH A BALLISTIC PISTON COMPRESSOR 



ENGINEERING AND INDUSTRIAL EXPERIMENT STATION 

College of Engineering 
University of Florida 
Gainesville 



SPECIAL REPORT 


THERMODYNAMIC PROPERTIES OF UF g 
MEASURED WITH A BALLISTIC PISTON COMPRESSOR 


By 

David E. Sterritt 
University of Florida 

George T. Lalos 
Naval Ordnance Laboratory 
White Oak, Maryland 

Richard T. Schneider 
University of Florida 


Work supported by 

National Aeronautics and Space Administration 
NGL 10-005-089 


March 1973 


Department of Nuclear Engineering Sciences 
University of Florida 



This work was supported by the National Aeronautics and Space 
Administration, N6L 10-005-089. 

The authors wish to thank Dr. Karlheinz Thom (US AEC/NASA, Nuclear 
Systems Office) for his help and scientific contributions to this work 
during his three month stay at the University of Florida. 

Thanks also go to Dr. Makoto Takeo, University of Oregon for 
supplying his simulation program and for helpful counseling in its use. 

Credit for the design and testing of the pulse shaping network 
goes to Mr. Gordon Hammond of the US Naval Ordnance Laboratory, Silver 
Spring, Maryland. 

The helpful advice of Dr. C. D. Kylstra, University of Florida, 
in the use of his optimization program is also acknowledged. 

We also would like to thank Dr. Richard D. Dresdner, University of 
Florida for his advice and help regarding the physics and chemistry of 



ABSTRACT 


From experiments performed with a ballistic piston compressor, 
certain thermodynamic properties of uranium hexafluoride are investi- 
gated. Difficulties presented by the nonideal processes encountered 
in ballistic compressors are discussed and a computer code BCCC (Ball is 
tic Compressor Computer Code) is developed to analyze the experimental 
data. The BCCC unfolds the thermodynamic properties of uranium hexa- 
fluoride from the helium-uranium hexafluoride mixture used as the test 
gas in the ballistic compressor. The thermodynamic properties deduced 
include the specific heat at constant volume, the ratio of specific 
heats for UFg, and the viscous coupling constant of helium-uranium 
hexafluoride mixtures. 
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I. INTRODUCTION 


Gas-driven pistons in closed-end tubes can be used to study the 
properties of gases. Fast response detection systems for determining 
the volume, pressure, and temperature of the test gas have made this 
area of investigation into the properties of gases most rewarding. The 
U.S. Naval Ordnance Laboratory pioneered research into this field in 
the early sixties [1]. Successful numerical analysis of the phenomena 
observed in ballistic compressors was achieved by Takeo [2] in 1965 with 
the development of a computer program which calculates gas leakage 
past the piston and heat losses from the test gas to the walls of the 
compressor. Measurement of the time lag between the maximum pressure 
and the maximum light intensity at the University of Oregon [2] supplied 
new evidence for explaining the phenomena of gas leakage, past the 
piston during compression. 

Using the initial work on numerical analysis done at the University 
of Oregon, the University of Florida has developed the BCCC which ex- 
tracts the thermodynamic parameters of the test gas from experimental 
pressures, temperature, and volumes, each measured as a function of time. 

The properties of uranium hexafluoride have been extensively in- 
vestigated at temperatures below 500° K. The specific heat ratios, the 
specific heat at constant volume, the vapor pressures, entropy, and 
enthalpy have all been investigated and tabulated but only at low tem- 
peratures. 

The specific heat ratios has been investigated by measuring the 
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velocity of sound and by measuring the ratio (Cp/C y ) at 50°C. This 
value is 1.063 ± 0.5%. 

Table 1 gives experimental measurements of the specific heat, C v , 
of UFg vapor. Measurements of this parameter have, by and large, been 
done by spectroscopic techniques. The result of these investigations 
have resulted in values for the specific heat at constant pressure 
which are 10% lower than other experimental measurements. 

Besides the obvious value of ballistic piston compressors in 
determining thermodynamic properties of gases at high pressures and 
temperatures, they can provide valuable information in other applied 
technology fields. The analysis of gases proposed for use in the gas 
core nuclear reactor is an excellent example of a low cost experiment 
to simulate conditions and analyze phenomena expected in large scale 
prototypes. In particular, the lifetimes of various species of gas 
molecules including that of uranium hexafluoride will be of extreme 
interest in the gas core reactor. 

Another application of ballistic piston compressors is development 
of the nuclear piston engine. By simulating one stroke of the nuclear 
piston engine, a ballistic piston compressor can provide valuable infor- 
mation needed before a critical nuclear piston engine is constructed. 

Also, by using a nuclear reactor in conjunction with a ballistic piston 
compressor, the interaction of fission fragments and plasma may be studied 
Since the plasmas obtainable with ballistic piston compressors are hot, 

dense, and long lived, they could provide a valuable research tool into 

) 

those areas where plasma lifetime has previously been a major restriction. 

Since the conditions of the test gas in a ballistic compressor 
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change by several orders of magnitude during a typical experiment, inves- 
tigations of gases at extreme conditions are easily feasible from initial 
conditions at standard temperatures and pressures. Maximum temperatures 
and pressures above five thousand degrees Kelvin and five thousand 
atmospheres are easily achieved. Thus, ballistic piston compressors 
are ideal devices with which to study gases under extreme conditions of 
temperature and pressure. 

Exact determination of the equation of state of the gas under 
investigation would necessitate continuous measurement of the temperature, 
pressure, and density of that gas. However, discrete measurement of 
these properties at even a selected part of the piston stroke provides 
enough information for the BCCC to analyze a substantial part of the 
equation of state. Thus, accurate information on the equation of state 
of a gas under investigation can be derived from measurements of the 
temperature, pressure, and density over part of the piston stroke. The 
best part of the stroke over which to obtain data is the end of the 
first compression stroke and the beginning of the first expansion stroke, 
thus including the maximum pressure and temperature, and the minimum 
volume. The measurements of these maxima and minima and their time 
relationships, in themselves, provide a wealth of information for analysis 

The experimental and analytical techniques developed in this report 
determine the thermodynamic parameters of the gas under consideration at 
high temperatures and pressures. To accurately describe the gas under 
consideration the Van der Waals' equation of state is used. One objective 
of these studies is the determination of the Van der Waals' constants 
a and b as functions of temperature. The ratio of the specific heats (i.e 
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constant volume and constant pressure) completes the thermodynamic para- 
meters needed for gas description. Since the two specific heats are 
related, only one need to be analyzed by the system. The specific heat 
at constant volume has been selected. 

Gas leakage past the piston, as shown later, is an important 
mechanism in ballistic piston compressors. In particular, the viscosity 
of the gas under test determines the leak rate. Since the leaking 
does not occur under conditions of streamline flow, the viscosity which 
determines gas leakage in a ballistic piston compressor will be a compo- 
site of the kinematic viscosity, dynamic viscosity, and associated turbu 
lent effects. Dissociation is negligible for the temperature and pres- 
sures under consideration [4]. These facts must be considered when 
analyzing the measured viscosity. 



II. EXPERIMENTAL APPARATUS 


A. Ballistic Piston Compressor 

The ballistic compressor used in the system is shown in Figure 1, 
and consists of four main parts, the reservoir, the piston release 
section, the tube, and the high pressure section. The compressor rests 
on five supports, the center three are steady rests, while the outer two 
consist of roller bearings. During assembly and alignment of the com- 
pressor all five supports are utilized. For firing, however, the three 
center supports are lowered and the compressor is then free to recoil on 
the two outer roller bearing supports without transmitting any force 
to the two tables. The two tables and the center support assembly are 
adjustable in height so that the compressor can be raised or lowered 
as necessary to align the radiation viewing window with the optical 
system. 

The reservoir, shown in Figure 2, was designed for a maximum 
operating pressure of 136 atmospheres and was tested to 200 atmospheres. 
The testing and certification was made by a separate concern from the 
one that fabricated the reservoir. With this reservoir it is possible 
to generate a 10,000 atmosphere test gas pressure for any test gas with 
an initial pressure of up to two atmospheres. The volume of the reser- 
voil is about 60 liters, sufficiently large that the driver gas pres- 
sure does not drop excessively as the piston travels down the barrel. 

The piston release section serves two functions: it fires the 

piston and in addition makes possible the sealing off of the reservoir 
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Figure 2. Ballistic Piston Compressor (Reservoir End) 
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gas from the tube in such a manner that the tube and high pressure sec- 
tion can be disassembled safely while the reservoir is fully charged 
with driver gas pressure. This feature greatly facilitates the experi- 
mental work, since it makes possible the exchange of pistons, windows, 
and pressure gauges, and allows inspection of the barrel bore and high 
pressure section, without time consuming and expensive dumping of the 
reservoir pressure. 

Figure 2 shows the piston release section in position between 
the tube and reservoir. It consists of a cylinder closed on one end with 
a cap to which is welded a length of pipe. This pipe goes through the wall 
of the reservoir to the gas supply system. Four ports are located in 
the walls in the front part of the cylinder. The two ends of the plunger 
when located in the front part of the cylinder bridge the four ports. 

The plunger necks down from a 7.61 cm diameter on the reservoir side to 
a 4.44 cm diameter on the barrel side. The large diameter part of the 
plunger contains two neoprene 0-rings while the small diameter part con- 
tains one. In moving to the rear position the plunger fires the piston. 

A safety cap is screwed onto the barrel end of the piston release section 
whenever the barrel or high pressure section is disassembled while the re- 
servoir is charged with driver gas pressure. An accidental loss of gas 
pressure behind the plunger, releasing the driver gas, would thus not 
endanger personnel or equipment. 

The barrel is 3.89 meters long and has a 0.09 mm bore. The material 
is A.l.S.l 4340 steel with a yield strength of 12,650 Kg/cm^; the hardness 
is Rockwell C 42-45. The 50 mm bore section is 17.78 cm long; the wall 
thickness is 2.57 cm. A gas inlet is located in the wall to facilitate 
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movement of the piston while checking the apparatus. This section was 
designed for pressures up to 5,000 atmospheres. Since the high tempera- 
tures which are generated in the high pressure section shown in Figure 
3 cause a slight vaporization of the inner walls with the resulting 
emission of unwanted impurity radiation, the inner bore was plated with 
chromium. 

Figure 4 shows the piston and cup seals in detail. To minimize 
gas leakage around the piston, two cup seals are employed. The gases 
before and behind the piston expand these seals against the bore thus 
minimizing gas leakage around the piston. The piston body is made of 
steel. Two phosphor bronze bearings provide the surface upon which the 
piston moves. The head of the piston is made of molybdenum to prevent 
ablation by the hot gas. 

The gas handling system is shown in Figure 5. Most of this system 
is housed in a mobile control console. This console controls the test 
gas and reservoir fill pressures. It also contains the vacuum system 
and helium purification system. Firing and post-firing procedures are 
implemented from this console. A relationahip between the reservoir 
pressure and the maximum pressure of the test gas can be easily derived. 
This is accomplished by equating the work done by the reservoir gas done 
by the piston to the work done by the piston on the test gas. The as- 
sumptions for this derivation are 

1) ideal test gas 

2) no gas leakage across piston 

3) adiabatic process 

4) constant reservoir pressure 
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Figure 3. Ballistic Piston Compressor (High Pressure End) 
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Figure 4. Piston and Cup Seals 
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Figure 5, Gas Handling System 
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This gives the following equation 


max 


- (-T 1 ) 


1/Y 


res 


<W /Y - i 


0 ) 


This equation shows that large values of maximum test gas pressures 
can be obtained for relatively modest reservoir pressures. 


B. Diagnostics 

1. Measurement of Test Gas Pressure 

The data acquisition system is shown in Figure 6. Possible 
optical misalignment caused by compressor recoil is monitored by means 
of a slidewire resistor, sensing the position of the high pressure test 
section. 

The pressure of the test gas is measured continuously with a cali- 
brated Kistler pressure transducer. A charge amplifier provides bias 
to the transducer. The oscilloscope trace of the output of the pressure 
transducer is photographically registered. The charge amplifier also 
provides a calibration signal for the range of interest, as well as a 
zero pressure reference. Both of these traces are photographed before 
each experiment. Figure 7 shows the results of a typical pressure mea- 
surement showing the pressure signal, the 5000 psi calibration trace, and 
the zero baseline trace. 

2. Measurement of UFg Pressure 

The pressure of the uranium hexafluoride is measured with a 
Wallace and Tiernan model FA160 absolute pressure gauge. Initially the 
test section is evacuated and then filled to the appropriate pressure of 
uranium hexafluoride. 
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Figure 6, 

DATA ACQUISITION SYSTEM 









Figure 7, Baseline , Calibration/ and Pressure Traces 





17 


The gas line from the uranium bottle to the compressor barrel con- 
sists of twenty inches of quarter inch copper tubing followed by four 
feet of one inch flexible tubing leading to the barrel. There is a 
quarter inch port (or "touch-hole") drilled in the barrel, just ahead 
of the initial piston position, to which the flexible tubing is connected. 

The tap to the absolute pressure gauge consists of three feet of 
quarter inch copper tubing which is connected to the flexible tubing. 

After each filling of the barrel and test section, but before each 
reading of the uranium hexafluoride pressure, at least five minutes were 
allowed for equilibrium to be established. 

After filling with uranium hexafluoride, the test section is 
filled with helium to a total pressure of one atmosphere. The total 
pressure is determined with a differential pressure gauge. 

3. Volume Measurement 

To measure the volume of the test gas mixture near maximum pres- 
sure, the high pressure test section is instrumented with a magnetic 
pickup, which senses the position of the piston during compression [5]. 

The body of the piston has circumferential grooves at measured 
distances. The voltage signal induced in the magnetic pickup by these 
grooves is applied to a pulse shaping network which produces a square 
wave for each groove which passes the magnetic pickup. The return of 
this square wave to baseline coincides with the middle of the associated 
"tooth". Figure 8 shows the voltage signal from the magnetic pickup 
and the associated square waves. The square wave generator also provides 
a radio frequency triggering signal beginning with the first groove. 

This signal is used to trigger the data acquisition system. Figure 9 
illustrates the geometric configurations of the piston and magnetic pickup. 
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Figure 8, Trace of Magnetic Pickup and Amplifier Output 
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Figure 9. 

DIAGRAM ILLUSTRATING MAGNETIC PICKUP 
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4. Temperature Measurement 

The temperature of the test gas is measured by the brightness- 
emissivity method. In this technique the emission and absorption of 
the gases are measured and the ratio of these values is used to determine 
the temperature of the test gas. 

Ideally the emission and absorption should be measured simultan- 
eously. However, they can be measured alternately if conditions in 
the test gas do not change appreciably between measurements. The latter 
method is used as shown in Figure 10. 

Because of the relatively low temperature of the test gas, the 
helium is not sufficiently excited to allow the measurement of tempera- 
ture to be made on a spectrum line of helium. Therefore a small amount 
of NaCl is introduced into the high pressure test section. The tempera- 
ture is measured using the 589nm sodium line. 

The He-Na collision cross section is sufficiently large that ap- 
proximate equilibrium is set up between the temperatures of the helium 
and sodium atoms despite the low electron density. 

The temperature of the test gas is given by: 

V V 3 

1/T = 1/T l + k/hv 1 n ( 1 + (2) 

where the following symbols have been used: 

k the Boltzmann constant 

h the Planck constant 

v the central frequency of the sodium line 

T-j the brightness temperature of the standard lamp 

V, the photomultiplier output voltage when the 

standard lamp is cut off by the rotating disk, 
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and hence corresponds to the intensity of the 
radiating sodium line 

Vp the voltage of the photomultiplier when only 
the standard lamp is emitting (i.e., with a 
cold plasma) 

Vg the voltage of the photomultiplier with 
standard lamp and emitting test gas 

A tungsten standard lamp is used. The lamp temperature is adjusted 
to be close to the temperature expected in the test gas. Figure 10 
shows the trace of the raw data of a typical temperature measurement. 

The validity of the above method hinges on the assumption of 
thermodynamic equilibrium of the system. Since the velocity of the 
piston is on the order of a tenth of the speed of sound, it is reasonable 
to assume that the gas is in thermodynamic equilibrium. 

Only the center portion of the pressure broadened sodium doublet 
was used for these measurements. This ensures that a high signal to 
background (continuum) ratio was achieved. Figure 11 shows the results 
of the temperature measurements using the brightness emissivity method 
compared to the temperatures computed by the BCCC. The measured values 
are higher than the computed values by about 25%, which is consistant 
with previously reported results [6], 



Figure 10, Raw Data of Brightness Emissivity Temperature Measurement 
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Figure 11, Experimental and Calculated Temperatures Vs. Time 
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III. THEORETICAL REVIEW 


A review of the theory employed in obtaining the numerical algori- 
thms used in the BCCC (Ballistic Compressor Computer Code) is presented. 
The Van der Waals 1 equation of state and the expression for the specific 
heat ratios derived from that equation of state provide introduction 
for their subsequent use in slightly modified form. 

Non-ideal processes are examined and reference is made to analyti- 
cal evidence suggesting that such processes cannot be neglected. The 
rate equations for gas leakage and heat loss complete the analytic 
expressions of primary interest. 

The numerical equivalent of the heat equation for the cylindri- 
cal bore introduces similar equations in the following chapter. The 
expression for the work done on the piston by the leaking gas and an 
empirical formula for the viscosity complete this. chapter. 

A. Equation of State 

The Van der Waals' equation of state is: 

(p + a/v 2 )(v - b) = RT (3) 

where a and b are the Van der Waals' constants. In general a and b 
are functions of temperature, and are represented by a second degree 

(4) 

(5) 


polynomial in temperature: 

2 

a = a-| + a 2 T + a^T 
b = b 1 + b 2 T + b 3 T 2 
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Theoretical temperature independent values may be obtained for 
a and b from "real gas" theory [6]: 


a - 3p c v c 2 

(6) 

b = v c /3 . 

(7) 


p and v are respectively, the critical pressure and the critical molar 
0 c 

volume. Table 2 lists the Van der Waal s' constants for helium; the 
values for uranium hexafluoride are calculated from its triple point. 


B. Specific Heat Ratios 

One of the principal objectives of this investigation is the 
determination of the ratio of specific heats of uranium hexafluoride. 
The specific heat ratio or "gamma" as it is sometimes referred to given 
by: 


Y = 


3 

R pv +av 
C v pv^+av-2a(v-b) 


(8) 


where C y is the specific heat of the gas at constant volume. In general 
this parameter is a function of temperature. For the uranium hexafluoride, 
a second degree polynomial representation in temperature is made of C y . 

The specific heat of helium is assumed to be constant over the range of 
conditions to be investigated. 

C. Non-Adiabaticity of the Compression Process 

Ideally, the ratio of specific heats would be measured in a poly- 
tropic process. That is, the pressure, temperature or volume of the 
gas of interest would be measured as the gas is compressed without loss 
of heat or mass. These conditions are approached in the rapid compression 
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of a gas by a piston. However, it has been shown [2] that while a gas 
is compressed by a ballistic piston, heat is lost to the barrel walls 
and gas is lost past the piston. Additionally, it would be expected 
in the ideal compression of a gas, that the maximum pressure and minimum 
volume would occur simultaneously. Results show they do not. As 
shown in Figure 12, the minimum volume may occur after the maximum 
pressure by as much as a millisecond. This phenomenon is caused by 
leakage of the test gas past the piston. It has been estimated that 
the leakage rate is on the order of ten moles of test gas per second 
at the point of maximum compression [2]. Clearly, this represents a 
significant departure from the ideal compression of a known amount of 
gas. 

The amount of heat lost from the test gas during the compression 
process is also significant. This non-adiabaticity may be examined by 
comparing the heat conducted through the wall to the heat necessary 
to raise the temperature of the gas during the compression stroke [2]. 

To assume adiabatic compression, the following inequality must 

hold: 

<«lost « HC v AT < 9 > 

The right hand side of this equation is the heat necessary to raise the 
temperature of the test gas by AT. After solving for the motion of the 
piston, the above inequality gives the relation [2]: 

o x n- x i 
2 ttK(X 0 -X) 2 H-Xp-r" 1 

M C y ln(l+a/y 0 ) 

« (y-1)v®Y 75T (X/X 0 ) 1/2 Xa n (X/X 0 ) n 


00 ) 


Figure 12 . 

TIME COMPARISON OF MINIMUM VOLUME a MAXIMUM PRESSURE 
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where the following symbols have been used: 

K the thermal conductivity of the barrel 
wall (steel) 

Xq the initial length of the test gas section 

X length of tube occupied by the test gas 

y specific heat ratio (for this equation 

the test gas is the same as the reservoir 
gas) 

M mass of the test gas 

C y specific heat of the test or reservoir gas 

a cylinder wall thickness 

r Q radius of inner barrel surface 

P ro initial pressure of the reservoir gas 
m mass of the piston 

a coefficients of the power series resulting 
from the solution for the velocity of the 
piston 

A cross sectional area of the piston 

The non-adiabaticity of the compression process is illustrated by con- 
sidering the end of the first piston stroke where the piston's velocity 
is nearly zero. This is where the test gas temperature exceeds room 
temperature by the greatest amount. It has been shown [2] that in this 
region the right hand side has a value of nearly zero while the left hand 

side has a value of about five. Thus, the above inequality is not 

satisfied and the process cannot be assumed to be adiabatic. 

The other significant departure from an ideal process is the mixing 
of the test and reservoir gases. This mixing is the result of leakage 

due to the pressure difference between both ends of the piston as well 
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as to the motion of the piston. The leaking rate is related to the 
viscosity of the leaking gas by the equation [2]: 

dn/dt = if r o ( V r p )[ 87;S- (r o- , 'p> 2 + <"> 

where the following symbols have been used: 
dn/dt leaking rate 

p density of the leaking gas 

m a molecular weight of the leaking gas 

r Q inner radius of the barrel 

r p radius of the piston 

n viscosity of the leaking gas 

V piston velocity 

3P‘/3X the pressure gradient along the piston 


The quantity (r Q -r ) is also known as the piston gap and, as the term 
implies is the clearance separating the piston and tube wall. 

Heat loss from the test gas to the bore wall must be calculated 
numerically and the energy thus lost must be subtracted from the inter 
nal energy of the test gas mixture. In cylindrical coordinates, the 
heat equation is: 


3T _ #3 • , 1 3T \ 

W ~ ? + ) 

3r 4 r 3r 


( 12 ) 


where a is the thermal diffusivity of the steel barrel. For numerical 

purposes, the bore thickness is divided up into nodes, and the tempera 

ture at each node is calculated. Let the subscript i denote the ith 

iteration in time, and the subscript j denote the node at radius 

[R. + Arl: 

L inner J 
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T i+1 »J~ T i >3 = JiJ+r CT 1;J+ T 1,J-l i , 1 f T i J+l" T i ,3-1 

At 1 Ar 2 } R inner + ^ 2Ar 

(13) 


where R. is the inner radius of the bore. The new temperature at 

inner r 

time step i+1 and node j becomes: 


T i+ i 4 = T. i+1 (-^2 + 1 )<*At + T i 

i i.J I A / 2(R+jAr)Ar 1,0 Ar 2 


+ T, 4 ,( 


1 


1 


Ar 2 2(R+iAr)Ar 


)aAt 


(14) 


where At is the width of the current time step. 

The work done by the leaking gas on the piston is given by: 

2RL Vn 

E = E_ttAX - (P -P)6nrAX (15) 

G r 

where L p is the length of the piston and G is the piston gap. 

The above expressions are implemented numerically to simulate a 
ballistic piston compressor. By taking short time steps (i.e., one 
microsecond), the motion of the piston may be calculated and the changing 
thermodynamic conditions monitored. That is, as the piston "steps" 
down the tube numerically, the changing velocity, temperature, and pres- 
sure, etc., may be accurately followed. 

Since several parameters needed at each time step change with 
temperature, etc., they must be recalculated. Included is the specific 
heat at constant volume of the UFg and the Van der Waal s' constants. 
Others must be recalculated as the ratio of UFg/He changes due to gas 
leakage. In fact, the characteristics of the test gas mixture may change 
dramatically due to helium leaking into the test gas at the beginning of 
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the compression stroke. Therefore, the thermodynamic parameters of the 
test mixture must be re-evaluated at each time step. 

An empirical formula has been used to describe the change in the 
viscosities of the gas as a function of temperature: 

n = AT B ' (16) 

A and B are also assumed to be linear functions in temperature: 

B-i +B 9 T 

n = (a-, + A 2 T)T 1 6 (17) 

Using this equation, the viscosity is also recomputed each time step. 



IV. ANALYTICAL TECHNIQUE 


Implementation of the theory is achieved by means of a numerical 
scheme to simulate the motion of the piston in a ballistic piston com- 
pressor. In this scheme, the movement of the piston is divided into 
steps and at every step, the nonlinearities are approximated and intro- 
duced into the equations of state. The final result of the simulation 
is a tabulation of the variables of interest as a function of time. 

In general the outputted variables include the length of bore occupied 
by the test gas, the ratio of specific heats of the test gas, the pres- 
sure, temperature, number of moles of the test gas mixture, the pres- 
sure of the reservoir, and the acceleration of the piston. 

The numerical and experimental data are then aligned in time and 
errors are assigned to the numerical data. By taking the partial deri- 
vatives of the errors with respect to those thermodynamic parameters 
being optimized a quantitative measure is obtained of how to adjust 
those thermodynamic parameters to reduce the error. When the error 
becomes insignificant, those thermodynamic parameters which accurately 
describe the gas have been found. 

A. Piston Motion Simulation 

The program developed to analyze the piston motion in a ballistic 
piston compressor models the first compression and expansion stroke of 
the compressor under analysis. As explained previously, there are 
several nonlinear processes which must be included in any accurate analysis. 
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The first requirement for the program then, is that its time step must 
be small enough to give adequate treatment to the nonlinearities such 
as gas leakage and heat loss. 

Each time step calculates piston movement, heat loss to the cylinder 
walls, gas leakage past the piston, and work done on the piston. In 
addition, since a test gas mixture of two species is to be analyzed, and 
since the proportions of the mixture may change with leakage, each 
time step must recalculate thermodynamic parameters of the test gas 
mixture. Figure 13 gives the flow chart of the simulation part of the 
BCCC. 

Each time step is begun by calculating piston displacement from 
the velocity and the width of that time step. 

AX = VAt (18) 


where X is the length of bore occupied by the test gas and V is the 
velocity of the piston. New temperatures and pressures are calculated 
on the basis of this change in volume: 


_ Y 2a(V-b) 

p i = ( p i_i + -%■)(! - 0± tttt) [(PV 2 +a)(V±AV)-2a(V-b)]C 

1 1 1 r V D v (V±AV) 2 

(19) 

_ Y _ £a.( V-b) 

(V±AV-b)(P._ 1 + ^-)(1.0± ^) [(PV 2 +a)(V±AV)-2a(V-b)]C v 

T. = - 

R-n 

( 20 ) 

Pj new pressure of current time step 
P i_l pressure of previous time step 
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a Van der Walls' constant multiplied by the 
number of moles of gas squared 

b Van der Waal s' constant multiplied by the 
number of moles under consideration 

V reservoir or test volume under consideration 

Y specific heat ratios of gas under consideration 

C y specific heat of gas under consideration 

R ideal gas constant 

n number of moles of test or reservoir gas 

under consideration 

± + for reservoir, - for test gas 

Gas leakage past the piston is proportional to the pressure dif- 
ference between the ends of the piston. Also the friction work of the 
leaking gas done on the piston is calculated at this time as well as 
the internal energy the leaking gas carries with it. 

The sign of the pressure difference is tested to determine direction 
of leakage and then the number of moles leaked is calculated from: 

Am = (- ^p" P ^ 2 + v) £ AT - G (21) 

oP^ri v 

The energy which the leaking gas carries with it is calculated by: 

E £ = (TC y + Pv)Am (22) 

and the energy lost from the leaking chamber due to friction work is given by 

2RL Vn 

E fw = [(Pr-P)rG - — ]ttAx (23) 

These are used later in the time step to calculate pressure cor- 
rections in the reservoir and the test gas mixture. 
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With these pressures, temperatures, volumes, and gas leakage a 
new evaluation of the thermodynamic parameters of the test and reser- 
voir gases is performed. This is necessary since the initial leakage of 
the reservoir gas into the test chamber can significantly alter the pro- 
perties of the test gas. In the present case, the test gas is made up 
of helium and uranium hexafluoride. A numerical composite of the pro- 
perties of these two gases is made at each time step. 

Approximating the uranium hexafluoride-helium mixture by a single 
species test gas is done by weighing the thermodynamic parameters of 
each gas by the respective number density of that gas. By this method, 
the initial conditions of the single species is obtained, and those 
thermodynamic parameters which change during the compression or expan- 
sion of the test gas mixture are updated. 

The number of moles in the single species approximation is equal 
to the sum of the moles of uranium hexafluoride and helium present. 

Those parameters which are updated at each time step, during com- 
pression or expansion include the Van der Waal s' constants, the specific 
heat at constant volume, and the viscosity. Also the number of moles 
of uranium hexafluoride and helium is updated at each iteration and 
the running totals are used in the update. When gas leaks into the 
test section from the reservoir side it is assumed that only helium 
has leaked into the test gas section and only the running total helium 
is updated. When gas leaks out of the test section, it is assumed that 
both helium and uranium hexafluoride have leaked out in proportion to 
their molar fractions present at the last time step. The relation used 
is: 
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n UFg = n UF 6 ^°^ d ^ ^ + n (lost) /n T(oldp 
n He = n \tel°' d W + n lost /n T(old) ) 


(24) 

(25) 


where n(old) indicates values from the previous time step, and n (-| ost ) 
is the number of moles leaked since the last time step. Since the opti- 
mization part of the BCCC assumes continuous parameters, the viscosity 
of the single species approximation of the test gas mixture is evaluated 
next. The viscosity coefficients "A" in equation 16 for the single 
species approximation are calculated from: 


A 1S ~ ^ n He /n T^A-| He + C l( n He n UF 6 /n T 2 ) 
+ ( n uF 6 /n T^ A lUF 6 


(26) 


A 2S “ ^He /n T^A 2He + C 2( n He n UF 6 /n T^ 


+ (n UF /n T ) A 

b 


2UF, 


(27) 


where the subscript "S" implies single species, n Hg , nyp , and n^ are 

6 

respectively the number densities of helium, uranium hexafluoride, and 
the total of the two [8]. The variation of the exponent, B for the 
single species approximation is assumed to be the same as that of helium. 
The final result for the single species viscosity is then: 


n$ = (A ls + A 2 S T)T 


( B lHe + B 2He T ) 


(28) 


The specific heat at constant volume of the uranium hexafluoride 
is approximated as a function of temperature by: 


c v (uf 6 ) = c, * c 2 t + c 3 T 2 


(29) 
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The specific heat at constant volume of helium is assumed to be constant 
in those regions of interest. Table 3 gives the specific heat coeffi- 
cients for helium and uranium hexafluoride. The specific heat at con- 
stant volume of the single species approximation is calculated by: 


'VS 


n He r " UF 6 r 

n T u VHe n T b UF 6 


(30) 


Because the Van der Waal s' constants of the helium and uranium 
hexafluoride cannot be combined in a simple relation to give coefficients 
of a single species, no updating is performed with respect to the mole 
fractions of helium and uranium hexafluoride on the Van der Waal s' 
constants. The Van der Waals 1 constants of the single species approxi- 
mation can initially be taken to be the same as those of helium, and 
then be optimized. Updating of the Van der Waals 1 constants is only 
done with respect to temperature. The relation used for the single 
species Van der Waals 1 constants are: 

a S = a lS + a 2S T + 3 3S t2 
b S = b lS + b 2S T + b 3S T 


As demonstrated previously, the ratio of specific heats (equation 
8) is a function of the specific heat at constant volume and the Van der 
Waals' constants. Since the Van der Waals' constants are not updated 
with respect to the mole fractions of the test gas, the ratio of speci- 
fic heats is to a degree uncoupled from the changing mole fractions of 
the test gas mixture. However, this degree is small since the specific 
heat at constant volume is updated (both temperature and composition) 
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Table 3 


Specific 

Heats at 

Constant Volume (C y ) of He 

and UFg 

Vapor 

UNITS 

(C v /R) 

(C V /R°K) 


(C V /R°K 2 ) 


C 1 

c 2 


c 3 

He 

0.125 

- 


- 

UF 6 [3] 

1.140 

8.21 x 10" 4 


- 


R— Ideal Gas Constant 
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and its effect on the ratio of specific heats is much greater. 

Heat loss to the bore walls and end plugs (one of which is the 
piston) is calculated each time step. The thicknesses of the cylin- 
der walls and end plugs are divided up into nodal points at which the 
temperature is evaluated. At most, a hundred nodes are used for either 
the cylinder walls or end plugs. The end plugs are assumed to be twice 
as thick as the cylinder walls. The formulae used are: 


T cj (t i ) = T c,j+l (t i-l )C ~2 + 


1 


Ar 2(R+(j-l ) Ar)Ar 


-]aAt 




cj 

+ T cj-l (t i-l )[ 7J “ 1 


Ar 2(R+(j-l ) Ar)Ar 


-]aAt 


T cl<V + H4rT sas (t,)/(l + H4r) 


W ’ T 


room 


(33) 

(34) 

(35) 


where the subscript j is the nodal index, t^ is the ith time step, and 
T c is the temperature of the cylinder wall. 

For the end plugs, the formulae becomes: 


aAt 


WVi )>- (2Ar)2 


+ T ,(t, ,)(1 - 2 

eJ 1 ' (24r) 2 

T el<‘i) ■ W + HArT ga s (t i )/(1 + HAr) 


(36) 

(37) 


The net heat flux into the test gas chamber is given by: 

2irrXk(T 2 ( t, )-T 2 (t 1 ) ) (T „,(t,)-T .(t,)) 

H- = ( ! £i_L_ + 2AK — ! “_L_) A t (38) 


Ar 


2Ar 
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The reservoir and test gas pressures are now adjusted to compen- 
sate for the heat flux out of the test gas, gas leakage, and work done 
on the piston: 


fw- 


P new ■ p + <>W { o » R/(v - b)C v 


Ynew ' P r + R /<V b r> C vr 


( 39 ) 

( 40 ) 


where the symbol (J^servoiMeakage^ has been used to differentiate 9 as 
leakage between the two chambers. 

The change in the velocity of the piston is now calculated from: 
(P r -P)AAX ♦ E fw 


Av = 


v 


( 41 ) 


The change in the velocity of the piston concludes evaluation of 
physical parameters for the time step. The final computation of the 
time step concerns the width: of the time step itself; whether it is to 
be increased, decreased, or remain the same. Two indicators are used 
in making this decision. They are the magnitude of the fractional change 
of velocity, i.e., | ( AV/V ) j and the fractional change in the pressure 
of the test gas mixture (P(t^) - P(t^_^ )/P(t^ ) ) . After consulting these 
indicators the time step is either increased by a factor of 1.2 or 
decreased by a factor of 0.8 or remains the same. 

Also at this point output of the time step is performed. Output 
variables of interest may include the time, pressures, temperatures, 
thermodynamic parameters, etc. 

Transfer is now made to the beginning of the next time step to 


repeat the process. 
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B. Numerical Error Analysis of Piston Motion Simulation 

The analysis of the data from the ballistic piston compressor con- 
sists of using the experimental data as target optimization by the BCCC. 

The BCCC takes as input initial guesses of the thermodynamic 
parameters, and the experimental results of pressure, volume, and tem- 
perature as functions of time. Computational use of the experimental 
data is made after each simulation of the motion of the piston and 
physical state of the gases in the ballistic piston compressor. The 
experimental data is used to assign error criteria to the numerical si- 
mulation. Since different parameters will be compared for error, frac- 
tional errors will be used. 

The origins in time of the numerical and experimental data are 
different. The origin of the experimental data coincides with the piston 
entering the high pressure test section. The origin of the numerical 
simulation coincides with the initiation of the firing sequence. The 
experimental and numerical data must be aligned in time before they can 
be compared. The offset between the two origins is approximated by 
considering the experimental and numerical minimum volumes and maximum 
pressures. The first offset is the difference in time between the ex- 
perimental and numerical maximum pressures. The second offset is the 
difference in time between the experimental and numerical minimum volumes. 
The final alignment is determined by taking the average of these two 
offsets. Fractional errors are then assigned to those calculated data 
points which correspond to experimental data points. If no calculated 
data point falls exactly on an experimental point, one is obtained by 
linear interpolation between the two adjacent calculated points. 
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The expression for the first alignment by maximum pressures is: 



(42) 


where Atp is the alignment offset between the experimental and the cal- 
culated times of maximum pressure, t and tp C are respectively the 
experimental and calculated times of maximum pressure. 

The expression for the second alignment by minimum volume is: 




(43) 


where the subscript v now implies minimum volume. 

The average alignment offset is then: 

At = (At p + At y )/2 (44) 

As explained previously, there is a time lag between the minimum 
volume and maximum pressure, the latter occurring first. 

The first error point is assigned to the difference in the calcu- 
lated and experimental time lags. 

e, ~vt ( ft -t ) - (t -t ))/(t -t ) (45) 

1 vc pc' ' ve pe''' ' ve pe' ' ' 

W is the weight assigned to this error point. 

The second error point is the fractional error of the calculated 
and experimental maximum pressures: 


e 2 = 2<% - p c >/p e < 46 > 

where P g and P c are the experimental and calculated maximum pressures, 
respectively. The third error point is the calculated and experimental 
minimum volumes: 
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2(V e - V c )/V e 


(47) 


where V g and V c are the experimental and calculated minimum volumes 
respectively. The weighting factor for these errors is two. 

Next fractional errors are assigned to the pressure and volume 
data points: 


e 3+n ^ P e(t n ) “ P c(t n )) /P e(t n ) j ,n=l,2,...,N 
e 3+N+j = ^ V e(t,) " V c(t,)^ /V e(t.) j 

J J J 


(48) 

(49) 


where there are N pressure points and as many volume points. 

Fractional errors are then assigned to the temperature points: 


'3+2N+m 


c (T e(t )‘ T c(t ) )/T e(tJ m=1 ' 2 H 


(50) 


~m' "'~m' m' 

The weighting factor for the errors assigned to the data points after 
the first three is one. 

A total error is also supplied the optimization program; 

3+2N+M o 


e = 


A 


(51) 


The optimization part of the BCCC, after several simulation runs, com- 
putes the numerical partial derivatives of the fractional error points 
to be optimized, with respect to the thermodynamic parameters. 

A zero derivative with respect to any parameter would imply that 
that parameter has no effect upon the simulation. In that case, no 
optimization could be performed on that parameter. While none of the 
derivatives will be zero, certain parameters will characteristically 
effect greater change in the simulation than others. In order of 
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decreasing sensitivity the parameters are 1) specific heat at constant 
volume, 2) viscosity, and 3) Van der Waals' constants. It is to be 
expected then, that the Van der Waals' constants will be the hardest 
to optimize, while the specific heat at constant volume is the easiest. 



V. DATA REDUCTION AND EXPERIMENTAL RESULTS 

Optimization theory pertinent to the BCCC is presented. The 
Gauss method of optimization with Lagrangian constraint [9,10] is 
used in the BCCC. To further introduce the mechanics of the BCCC, a 
very brief description of the input data required by the BCCC is presented. 

The BCCC is capable of compensating for conditions not explicitly 
incorporated, such as the irregular shape of the piston, non-streamline 
flow, or boundary layer effects. This compensation is achieved by 
calibration of the BCCC. 

A. Method of Optimization 

The process by which the thermodynamic parameters are determined 
is one of optimization. Initial values of these parameters are required 
and obtained by an initial estimate. The fractional errors of the cal- 
culated data points with respect to the experimental results are used 
to modify these parameters by an appropriate logic scheme. 

To be more useful , such a logic scheme for parameter optimization 
must have a limit on the hypersphere being searched, or in other words, 
a boundary must be set on the values the parameters may take. In the 
BCCC, these limits are in the form of a one dimensional radius about 
the initial value of each parameter P^. 

Pi ± AP. (51) 

where Pi is the initial value of parameter P.j . It follows then that 
the degree of the hypersphere being searched is equal to the number 
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of parameters being optimized. 

As explained in the previous chapter, the sum of the square of the 
errors is minimized, 

3+2N+M 0 

e = l e. 2 (52) 

i=l 

Taking the partial derivative of these errors with respect to each 
parameter being optimized gives 

0 = Ve = 2V e (P.)e T (P.+AP.) (53) 

where it is assumed that the derivative Ve remains essentially constant 
over the proposed step size. Expanding the right side gives 

0 = 2ve i (P i )[e T (P i ) + Ve T (P i )AP i ] (54) 

In order to prevent the coefficient of AP . , the step size, from 
getting too small forcing a larger AP^ (as yet of unspecified width) a 
constant is now added to impose a minimum size on the coefficient of 
AP i . 

0 = ve(P i )x e T (P.) + VeVe i T (P.)AP i - XlAP i (55) 

> 

where I is the unit matrix, and X is the minimum coefficient of AP... 
Solving for AP^ gives: 

AP i = -Ve(P i )e T (P i )[Ve(P i )Ve T (P i )-XT]" 1 (56) 

The method described above is the Gauss Method with Lagrangian 
Constraint. A standardized program exists at the University of Florida 
to implement this technique and is incorporated into the BCCC. 
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B. Organization of the Ballistic Compressor Computer Code 

Organization of the ballistic compressor computer code, in terms 
of input, consists of four tables. The first table provides constraints 
to the parameter optimization section, such as the maximum number of 
steps to be taken in the hypersphere. The second table contains all 
parameters which may be optimized. The third table specifies the dimen- 
sions of the search hypersphere. The fourth table contains constants 
which are not optimized, and the experimental data points. 

From this initial information, the search program calculates, by 
methods previously described, to the set of coordinates in the parameter 
hypersphere having the smallest error. The accuracy of this final 
coordinate is a function of the number of steps allowed the search 
program, the accuracy of the initial parameters, and the reliability 
of the experimental data. 

The accuracy of the numerical model must also be included in the 
error analysis. Table 4 gives the results of a numerical simulation 
with no gas leakage or heat loss. Also listed are the initial condi- 
tions of the simulation. The close correspondence in the ideal case 
satisfies the requirement that the BCCC in the limit of the ideal case 
must reproduce theoretical answers. 

C. Calibration of the BCCC 

Calibration of the BCCC makes numerical allowance for those pro- 
cesses which have not been considered explicitly in the BCCC. Included 
are the irregular shape of the piston, non-streamline flow, and boundary 
layer effects between the test gas and the barrel wall. 

The calibration of the BCCC is done by adjusting the numerical 
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Table 4 

Results of Numerical Simulation 
Without Gas Leakage or Heat Loss 


INITIAL CONDITIONS 

P(reservoir) = 15.3 atm 
P(Helium as test gas) = 1.00 atm 
P(UFg as test gas) = 0.00 atm 
T(Temperature of test gas) = 300°K 
v(Molar volume of test gas) = 24.62 (Jl/mole) 

Conditions at Maximum Compression 


Numerically Simulated 

% Error of Ideal 

P(max)...417 atm 

1.5 

T(max)...3319 °K 

1.5 

v(min)... .653 5,/mole 

1.5 
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value of the piston gap, the piston length, and the thermal coupling 
of the test gas to the barrel wall. Adjustment of the piston gap and 
length compensate for the turbulence and non-steaml ine flow of the 
leaking gas. Thus while the actual measured piston gap or length may 
serve as an initial guess, those values giving the most accurate results 
will in general be different. Calibration of the thermal coupling of 
the test gas to the barrel wall compensates for . turbulence of the test 
gas and/or boundary layer effects. These phenomena affect the amount 
of heat energy which is conducted to the wall. 

The calibration constants are roughly independent of the type of 
test gas used. Thus, once determined for specific physical conditions, 
they should not change for different test gases. The invariance of 
the calibration under change of test gas permits the separation of the 
thermodynamic and physical calibration parameters. By performing the 
calibration with a test gas of known thermodynamic parameters, the cali- 
bration constants are the only quantities under search and may be deter- 
mined by themselves. Once determined, the process is reversed, and the 
thermodynamic parameters are searched. 

D. Analysis of the Ratio of Specific Heat 

Two primary mechanisms govern change in the ratio of specific 
heats during an experiment. The mechanism of greatest effect in experi- 
ments is the mixing of the test and reservoir gases. Figure 14 shows 
the number of moles of gas in the test gas chamber during a typical 
experiment. 

The positive slope of the curve in this figure defines that period 
when reservoir gas leaks into the test gas chamber; negative slope 
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implies leakage of the test gas into the reservoir. Note that while 
reservoir gas is leaking into the test gas chamber for well over 80% 
of the time, the net leakage by the end of the compression stroke is out 
of the test gas chamber. This results from the very high test gas pres- 
sure relative to the reservoir pressure during the latter part of the 
compression stroke. 

Figure 15 shows an initial increase in the ratio of specific 
heats brought about by leakage of the helium driving gas into the test 
gas. This increase is terminated at the point of time when the direction 
of leakage reverses (i.e., when test gas pressure becomes equal to the 
drive pressure). 

When plotted against temperature as in Figure 16, the ratio shows 
the same sharp increase for the same reasons as in the case versus pres- 
sure. 

Figure 17 gives the ratio of specific heats as a function of time. 
After the initial increase due to the leaking helium, the ratio goes 
through a minimum at the point of highest temperature. This decrease 
in the magnitude of the ratio of specific heats is a result of the in- 
creased degrees of freedom excited by the high temperatures, reflected 
in a higher value of the specific heat at constant volume. 

Since the ratio of specific heats of helium is still constant 
at the conditions under investigation, the exhibited minimum is due to 
the uranium hexafluoride. Table 5 lists the experimental values of 
specific heat (C y ) and the viscous coupling constant (C-j). 

E. Analysis of the Viscous Coupling of Uranium Hexafluoride and Helium 

The viscosity of uranium hexafluoride is known from 0°C to 263°C [3]. 




Figure 16. Specific Heat Ratio Vs, Temperature 
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Figure 17. Specific Heat Ratio Vs. Time 
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Table 5 


Experimental Values of Specific Heat (C y ) and 
Viscous Coupling Constant (C-j) 


C v< UF 6> ‘ C vl + C v2 T + C v3 l2 
(C yl /R) = 13.89 R 

(C v2 /R°K) = R/°K 

(C v3 /R°K 2 ) = -2.1 0 -6 R/°K 2 

Viscous coupling constant of He-UFg mixtures: 
C-j = 70 ypoise 25 ypoise 


(see equation 26) 
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The measurements of these investigations, fitted to the empirical formulae 
previously described, are extrapolated to the temperatures needed in 
this investigation. The viscous coupling constant C-| is extracted by 
the BCCC from experimental data. The results thus far set C-| at 70 
micropoise with an estimated uncertainty of 25 micropoise (see Table 5). 

Figure 18 shows the viscosity of the helium uranium hexafluoride 
mixture as a function of time and in Figure 19 as a function of pressure. 

F. Results of Numerical Analysis of Experimental Data 

The departure from a non-polytropic process has been shown to be 
due to heat loss to the barrel walls from the test gas, and friction 
work done by the leaking gas. Figure 20 is the P-V diagram of the 
compression process. Straight lines with slope of 1.48 are included for 
comparison. The average value of the slope of the experimental curve 
is very close to 1.48, which is also the average value of the ratio of 
specific heats as shown in Figure 15. One may conclude from this, that 
the polytropic exponent of this compression process can be approximated 
very closely by the ratio of specific heats (C p /C v ) . 

The ratio of specific heats for pure UFg shown in Figure 21 as a function 
of temperature with pressure as parameter. The figure is computed from 
the measured values of the specific heat at constant volume and the Van 
der Waals' coefficients based on the values of the triple point of UFg. 

The dashed part of the lines in Figure 20 are not physically realized, 
since UFg is not a vapor at these conditions. 

The low temperature values (T<500°K) of the solid lines exceed the 
values given in Table 1 [3], by 6%, while for temperatures greater than 
500°K our values are approximately 10% higher. 




Figure 18. Viscosity of the Test Gas Mixture 
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Figure 19, Viscosity of the Test Gas Mixture as a Function of Pressure 


Figure 20, P-V Diagram of the Compression Process 
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VI. CONCLUSIONS 


A computer program (BCCC), capable of deriving certain thermody- 
namic properties of a test gas from measurements made in a ballistic 
piston compressor, was developed. This program was employed to deter- 
mine the C^ term (see equation 29) of the specific heat equation (C v ); 
the viscous coupling constant of the UFg (see equation 26), at experi- 
mental conditions beyond the range of previous experiments. During 
the course of these studies it was found that there is no significant 
change in the Van der Waal s' constants a and b of both He and UFg over 
the range of conditions of the experiments. 

The principal goal of this investigation, the measurement of 
specific heat ratio of pure UFg (see Figure 21) was accomplished by 
the determination of the temperature dependence of the specific heat 
at constant volume at high temperatures. Further analysis of the 
results showed, that for the He-UFg mixture ratios employed the poly- 
tropic exponent (- InP/lnV) can be approximated very closely by the 
specific heat ratio (C p /C v ). 

Finally, it is evident that the successful implementation of the 
BCCC as indicated by these studies has materially enhanced the usefulness 
a ballistic piston compressors for investigation of physical properties 
of gases of high temperatures and pressures. 
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APPENDIX A 


DETAILS ON THE COMPRESSION PROCESS 
Further details concerning the processes during compression can 
be obtained from Figures A-l , A-2 and A-3. 

Figure A-l show the rate of heat loss due to heat conduction to 
the barrel and end wall as a function of time. Figure A-2 shows the 
rate of heat loss as a function of pressure. Figure A-3 shows the 
rate of friction work done by the leaking gas on the piston as a function 
of time, according to equation 15. Positive friction work implies work 
done by the driving gas on the piston (during the time period while 
driving gas is leaking into the test gas). Negative friction work 
implies work done by the test gas on the piston (while leaking out of 
the high pressure section). 
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Figure A-2. Heat Losses as a Function of Pressure 
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